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We formulate a tractable calculation scheme for a finite temperature two-loop resummed effective 
potential in gauge theories without using a high-temperature expansion. We apply it to the Abelian- 
Higgs model and its extension to a minimal supersymmetric standard model-like model and study 
the thermal phase transition. It is shown that our scheme improves the previous results about 10% 
in the quantities relevant to the phase transition, and its impacts on bubble dynamics could be even 
more sizable. It still holds that the stop-stop-gluon sunset diagram enhances the strength of the 
first-order phase transition even without the high-temperature expansion. 
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Introduction. The finite-temperature (T) effective 
potential is a powerful tool for studying various phe- 
nomena in hot early Universe, for example, thermal 
phase transitions at quantum chromodynamics (QCD) 
and electroweak scales. It is known that a first-order 
electroweak phase transition (EWPT) is one of the nec- 
essary conditions for successful electroweak baryogenesis 
(EWBG) 0. A naive perturbation method at high T is 
known as unreliable. Generally, the perturbative expan- 
sion in a coupling constant is invalidated at high T due 
to the sizable T-dependent loop corrections. Therefore, 
such higher-order corrections should be resummed to ob- 
tain a sensible result. In the EWBG context, the EWPT 
up to two-loop order have been studied so far. However, 
since the calculation of the resummed two-loop effective 
potential is notoriously a formidable task, a high-T ex- 
pansion (HTE) in which a particle mass is less than tem- 
perature was exclusively used. Since the relevant param- 
eter region for the successful EWBG is vc/Tc > 1, where 
Tq denotes a critical temperature and vc is a Higgs vac- 
uum expectation value (VEV) at Tc, the validity of the 
HTE is far obvious. Indeed, it turns out that the HTE 
of the two-loop sunset diagrams is valid only for a small 
mass ratio to temperature, < 0.01 Q, rendering all ex- 
isting two-loop order results of the EWPT based on the 
HTE in literature less precise. 

In the minimal supersymmetric standard model 
(MSSM), it is known that vc/Tc at the one- loop order 
is not large enough to satisfy the sphaleron decoupling 
condition even at a bubble nucleation temperature Q. 
However, the deficit in vc/Tc can be compensated by 
a potentially sizable two-loop sunset diagram consist- 
ing of scalar top (stop) and gluon 0, Q. Therefore, it 
is of great importance that such a correction is quanti- 
fied on more stuff basis than the HTE. The viability of 
the MSSM baryogenesis depends not only on the current 
Large Hadron Collider data but potentially on the two- 
loop calculation of the EWPT beyond the HTE. For the 
former, refer to the recent papers 



In this Letter, we formulate a tractable calculation 
scheme for the two-loop effective potential and thermal 
resummation without relying on the HTE. In order to go 
beyond the HTE, in which T-dcpendent divergent terms 
are omitted, we must execute the renormalization with 
T-dcpcndent counter terms arising from the resumma- 
tion. As an example, we demonstrate our calculation 
scheme in the Abelian-Higgs (AH) model and apply it to 
a MSSM-like model for studying the thermal PT. 

Model. Let us begin by introducing the AH model 
and define the various notations used here. The La- 
grangian of the AH is 
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The scalar potential is given by 
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We parametrize the scalar field in terms of the VEV (v) 
and fluctuation fields 



$ = — f={v + h + ia), 
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where h is a physical state and a is a Nambu-Goldstonc 
(NG) boson which is eaten by the gauge boson. The 
field-dependent scalar and gauge boson masses are 
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where we work in the Landau gauge. Throughout the 
paper, we will adopt MS scheme for renormalization. 

Resummation method. We closely follow the resum- 
mation method presented in Ref. [7|. The basic proce- 
dure of this resummation is simply to add and subtract 
the T-dependent mass terms Aym 2 in the original La- 
grangian. 
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First, we focus on the mass term of h. The resummed 
Lagrangian and new counter terms in the d dimension 
are respectively given by 



V 



( -8v 2 - A T m 2 + ) h 2 + ^ (6) 



where e = 4 — d, Sv 2 and SX are the usual MS counter 
terms that are determined by the 1/e terms, and /Lt is a 
renormalization scale. The same procedure holds for the 
mass of a. With the resummed Lagrangian, the scalar 
propagators are 
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where m 2 h a (T) = m\ a + A T m 2 ha . 

The above resummation scheme has been applied to 
various models with and without spontaneous symmetry 
breaking. The proofs of renormalizability and NG theo- 
rem in O(N) 4> A theory are given in Ref. Q. 

For the gauge sector, similarly, we add and subtract 
a thermal correction in the original Lagrangian. In this 
case, however, we should treat longitudinal and trans- 
verse parts of the gauge boson propagator separately 
since thermal corrections to them arc different from each 
other. The resummed bare Lagrangian is given by 
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where £^ mBum 



is the unrcsummed bare Lagrangian and 
L^i/(p) and T fll ,(p) are the projection tensors which take 
the form of 
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in the rest frame of the thermal bath. Note that addi- 
tional T-dependent terms are nonlocal and noncovariant, 
but gauge invariance is still maintained. With Eq. ((5J, 
the resummed gauge boson propagator is cast into the 
form 

D^uip) = — — l —L flv (jp) + — — —rTp V (p) 
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FIG. 1. The diagrams in which the gauge boson is involved 
in the AH model. 
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FIG. 2. The IPI two-point functions of A„. 



where m 2 L T = m 2 A + A^to^ t and r is an arbitrary real 
parameter. Note that the noncovariant part in the last 
line of Eq. (fTTj) yields less ultra-violet divergent loop in- 
tegrals. Requiring g^ u (T^p) - rP llu {p)) = 0, one gets 
r = (d — 2)/(d — 1). We will see that this specific choice 
of r makes it easy to perform the two-loop calculations 
significantly. Let us denote the second and third lines 
in Eq. (HU by D^{p) = D c °J(p) + 5D^(p). With this 
gauge boson propagator, together with Eq. ([7]), we will 
compute the two-loop diagrams. Since there is no diffi- 
culty to obtain the resummed diagrams for purely scalar 
sector, we here concentrate on the diagrams involving 
the gauge boson as depicted in FIG. Q] Each diagram is 
expressed as 
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where the sum-integral symbol is defined by 
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with k° — ko = ioj n = 2niriT and f^"(g) is the onc- 
particlc-irrcduciblc (1PI) two-point function of A M eval- 
uated with the propagator D C °J . Now we expand above 
[i e Vn(v, T)'s in powers of SD^ as 
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A T m 2 L L{q) - A T m 2 T T{q) SD^(q) 



A h (k + q)SD^(k)SD^(q). 
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Note that the expanded terms to first order in derivative 
can be written in terms of T\ (q) with i = (a) — (d) which 
are represented in FIG.[5J For the covariant sector in the 
first line, one can show that all T-dependent divergences 
are cancelled out and renormalization is successfully car- 
ried out in the ordinary manner Q . Hence, we here prove 
finitencss of the noncovariant sector. 

In general, because of the gauge invariance the polar- 
ization tensor is written as 

IVte) = H (q 2 )P^(q) 

+ An L (g°, q)L^(q) + ATl T (q° , q)T^(q). 
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An LjT (g , q) in the infrared (IR) limit (q° = and q 
0) may take the form 
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Thus, in the IR limit the second and third lines of 
Eq. (fT8|) are reduced to 
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Therefore, we obtain the finite result. Note that the sub- 
traction terms Axm/^ T in the counter term play an im- 
portant role in this proof. It is easy to show the finiteness 
of this noncovariant part without the IR limit as well Q . 
Likewise, the last line of Eq. (|18[) is also finite It 
should be stressed that so far the HTE is not used. 

Before studying the PT, we look into some approxi- 
mations that are commonly used in the finite-T two-loop 
calculations. Let us define 
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FIG. 3. K—(a,a,0), 3A'(2a/3) and 3A HTE (2a/3) are plot- 
ted as a function of a. 



where ai — rrii/T and the momentum is defined in the 



Euclidean space, and so k 



k, q 2 



q, etc. 



K (aj, a,j , a*;) is defined by the function that contains 

two ub{E), where, for instance, riB(Ek) = l/(e Ek ' T — 1) 



with E, 
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We define K (a, a, a) 



K{a) [7|. The ellipses stand for the zero-T contributions 
and the terms including one ns{E), which can be ex- 
pressed in terms of the one-loop finite-T integrals. 

In literature, in addition to the HTE, the following 
approximation is frequently used [l(| 



H(mi,m2,m 3 ) -> H 



mi + m 2 + m-3 



(23) 



We investigate the validity of this mass-averaging approx- 
imation without the HTE. Here, we take a± — a 2 = a and 
as = 0. This particular choice corresponds to, for exam- 
ple, the stop-stop-gluon sunset diagram in the MSSM. 

The numerical comparison of K (a, a,0), 3A^(2a/3) 

and 3/v HTE (2a/3) are shown in FIG. H We can see 
that all of them would agree in the a — > limit. It is 
found that K—(a,a,0)/(3K(2a/3)) < 0.6 for a > 1. 
However, the corresponding value of each function is less 
than one. So the error of the mass-averaging approxima- 
tion can be small On the contrary, the use of the 
HTE gives a large deviation and yet negative value at 
a = 0(1) even though K(a) is positive-definite by def- 
inition. It is concluded that the HTE of K(a) is valid 
only for a < 0.01, which is consistent with the results in 
Ref. 0. 

Application to a MSSM-like model. Now we move on 
to discuss the thermal PT using the scheme formulated in 
this Letter. As a first step toward the complete analysis 
of the two-loop driven first-order EWPT scenario such 
as the MSSM, we consider an extended AH model in 
which additional U(l) gauge boson and complex scalar 
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FIG. 4. vc /To in the three cases are shown as a function of 
P3. The input parameters are given in the text. 
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FIG. 5. The shown is the effective potential at Tc- We 
take g-i = 1.2. The remaining parameters are the same as 
in FIG. H 



with the HTE. The one-loop calculation of vc /Tc is also 
shown by the dotted black curve. This figure shows that 
enhancement of vq/Tc due to the t-t-G^ diagram can 
still persist beyond the HTE. In this specific example, 
the use of the HTE leads to the underestimated vc/Tc- 

The height of the barrier between the two degenerate 
vacua in the effective potential (denoted by hb) is also 
relevant to dynamics of the first-order PT. The effec- 
tive potentials at Tc and g% = 1.2 in the three cases are 
plotted in FIG. [5] The color and line coordinates are 
the same as in FIG. |U We find that vq /7c*|i-ioo P = 
186.55/76.75 = 2.43, v c /T c \ 2 -i oop = 191.80/74.81 = 
2.56 and vc/T c | 2 _ loop H te = 204.94/81.32 = 2.52. The 
significant increase of Tc in the HTE case may be the 
consequence of the artificial negative contributions to the 
quadratic term in the scalar potential Q . It is also found 
that hb at the two-loop level is somewhat higher than 
that of the one-loop case, delaying the onset of the PT. 
However, we may get the overestimated result once the 
HTE is used. We observe that generally the larger 33 can 
bring the larger errors in vc, Tc and hb- 

Conclusion. We have formulated the tractable calcu- 
lation scheme for the finitc-T two-loop resummed effec- 
tive potential in gauge theories. Our analysis on the PT 
indicates that for the typical parameter set the use of 
the HTE in the sunset diagrams can lead to the over- 
estimated results by about 10% in Vc and Tc and by 
around 50% in hb- Since our calculation scheme is model- 
independent, it is applicable for any C/(l) models. Even 
in the non-Abelian case, we expect that the dominant 
effects relevant to the PT would be modeled by the sun- 
set diagrams in the AH or its extended models since a 
gluon-gluon-gluon type sunset diagram has little effect. 



are introduced. The added Lagrangian is 

ac = -\g^ + |ivT - i< + y 2 m 2 )\t\ 2 + ^T, 

(24) 

where = d ll G v - d v G^ and D^t = (<9 M - ig^G^t. 
As is done in the AH model, we carry out the thermal 
resummation in our scheme. To be specific, we keep only 
0{T 2 ) corrections in A T m 2 and tt^(T). 

As mentioned above, in the MSSM the stop-stop-gluon 
sunset diagram enhances vc/Tc- We scrutinize this ef- 
fect with and without the HTE. In what follows, by the 

HTE we mean that only K is evaluated using the 

HTE. 

In FIG. H vc/Tc is shown function of 53. We 
set v = 246, nth = 35, t?1q = 0, y = 1.0, e = 0.5, A = 
0.3 and Jx = l^e^^n = 150, where j E ^ 0.577. 
The red curve represents the two-loop calculation, and 
the blue dashed curve denotes the two-loop calculation 
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